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Abstract. The purpose of this paper is to show that, at least for Lagrangians of mechanical type, 
nonholonomic Euler-Lagrange equations for a nonholonomic linear constraint D may be viewed as 
non-constrained Euler-Lagrange equations but on a new (generally not Lie) algebroid structure on 
D. The proposed novel formalism allows us to treat in a unified way a variety of situations in 
nonholonomic mechanics and gives rise to a version of Neoether Theorem producing actual first 
integrals in case of symmetries. 
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1. Introduction 

I There are many approaches to geometric mechanics in the literature. We will work with a natural 

J> generalization of the framework for studying mechanical systems proposed by W. M. Tulczyjew [2l)l[3T)] 

t" — (see also [32] and references therein). In the simplest form, the phase dynamics of the system is 

understand as the lagrangian submanifold L of the symplcctic manifold (TT*M, dy^A/) equipped with 
tangent lift d-r^A/ of the canonical symplectic form ujm of T*M. Here, M represents the configuration 
manifold of the system and T is obtained from the lagrangian submanifold dL(M) C T*TM induced 
by the Lagrangian L : TM — > R via the canonical isomorphism em '■ T*TM — > TT*Af. In other 
words, the phase dynamics, as well as the Euler-Lagrange equations, are obtained in a simple way by 
means of the Tulczyjew differential Al = Sm ° dL : M — > TT*M . It is important to observe that 
both TT*M and T*TA/ arc double vector bundles over T* M and TM (see [T5] and references therein). 
The resulting submanifold L of TT*M is a particular case of modelling dynamical systems as implicit 
differential equations defined by differential inclusions (see [2"51 12"4"] 1. 

This framework admits an immediate generalization for more general morphisms e : T* E — * TE* of 
canonical double vector bundles associated with a vector bundle E (see [H1IE])) inducing the Tulczyjew 
differentials Al = e o dL : M — * TE*. This generalization includes as a particular case a theory 
of mechanical systems based on Lie algebroids, as proposed by A. Wcinstcin [36 and developed by 
many authors, however in a different geometrical setting (see, for instance, [6] [7] [19j [20l [25] [26]). The 
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motivation for study systems on Lie algebroids is that they often appear naturally as results of some 
reduction procedures. This is is a situation similar to the one known in the theory of Hamiltonian 
systems: reductions may lead from a symplectic to a Poisson structure. 

An additional challenge and one of the most fascinating topics in geometric mechanics is the study 
of constraints in this context. Of course, a general problem of putting constraints for the system in 
a variational setting involves constraints for velocities as well as constraints for virtual displacements, 
as was noticed already in |31j . In some cases, however, one assumes that the constraints can be 
determined from a constraint subset D of TM (or, of E in the algebroid case) by certain well-described 
procedures. The best known approaches of this type refer to the so called vakonomic and nonholonomic 
constraints. In the simplest situation, for D being a linear nonholonomic constraint, i.e. just a vector 
subbundle of TM (or, of E in the algebroid context), this procedure describes the nonholonomic Euler- 
Lagrange equations by means of the d'Alembert principle, having analogs also in the algebroid case 
[71 [T51 [57J [5]. We should stress that our nonholonomic constraints are linear in the broader sense, 
i.e. they are subbundles over submanifolds of the original base manifold. The nonholonomic Euler- 
Lagrange equations are commonly viewed as being not variational equations. In [8] it has been pointed 
out that it is not exactly the case, if we extend slightly our understanding of Variational Calculus. 

In this paper we continue studying nonholonomic constraints on algebroids and showing that, at 
least for Lagrangians of mechanical type, the nonholonomic Euler-Lagrange equations are just non- 
constrained Euler-Lagrange equations but for a special algebroid structure on the constraint subbundle 
(see also [IB])- This shows that mechanical systems based on general (not necessary Lie) algebroids 
appear naturally in the presence of nonholonomic constraints and gives a powerful geometrical tool 
when dealing with constrained systems. In particular, we get a version of Noether Theorem with true 
first integrals for nonholonomic systems. We do not get all possible algebroids on D applying our 
procedure. In particular, if the original structure was a Lie algebroid, then the new algebroid bracket 
is automatically skew- symmetric, so we deal with a quasi-Lie algebroid. One can the associate with the 
sequence of procedures, like reduction by symmetries and passing to a nonholonomic constraint, the 
sequence of the corresponding novel structures serving as appropriate geometrical tools in describing 
the systems: 

_ , reduction by symmetries nonholonomic constraint . T . , , . , 

TM > Lie algebroid > quasi-Lie algebroid . 

All this is of course closely related to the discovery of the role of the nonholonomic quasi-Poisson 
brackets [55J [33 [TSJ H] , this time not in the Hamilton but in the Lagrange picture. 

The paper is organized as follows. In the next section, we recall after [S] the basic ideas of developing 
mechanics on a general algebroid E, in particular, the Euler-Lagrange equations. Then, we construct 
in this setting an analog of the Tulczyjew differential for linear nonholonomic constraints, together 
with the corresponding nonholonomic Euler-Lagrange equations. In Section 4 we study reductions of 
a general algebroid to an algebroid on a nonholonomic constraint D (satisfying a natural admissibility 
condition) along a given projection. We discuss also the problem, what algebroids can be obtained in 
this way, if we start with Lie algebroids. 

Section 5 is the most important part of our paper. For Lagrangian functions L of mechanical 
type ('kinetic energy - potential'), we show that the non-constrained Euler-Lagrange equations on 
D, derived for the reduced Lagrangian I = Lm and for the reduced algebroid structure on D along 
the orthogonal projection associated with the kinetic energy, coincide with the nonholonomic Euler- 
Lagrange equations. Moreover, for this nonholonomic case, we can apply therefore the generalization 
of the Noether Theorem proved in [3] to obtain actual first integrals. Passing to the nonholonomic 
constraints does not requires therefore any change in our unified algebroid approach to mechanics: 
nonholonomic Euler-Lagrange equations are included in our framework. 

We end up with two well-known examples of nonholonomic constraints, the Chaplygin sleigh and 
the snakebord, to show how simply the corresponding equations of motions can be derived by means 
of our method. 
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2. Geometric Mechanics on general algebroids 

Let M be a smooth manifold and (x l ), i = l,...,n, be a coordinate system in M. Denote by 
tm ■ TM — > M the tangent vector bundle, with induced coordinates (x 1 ,^), and by ttm ■ T*M — > M 
the contangent bundle, with induced coordinates (x l ,Pj). 

Let t : E — > Af be a vector bundle and 7r : E* — > M its dual bundle. Taking a local basis {e a } 
of sections of E, then we have the corresponding local coordinates (x l ,y a ) on E, where y a (e) is the 
ath-coordinate of e € -E in the given basis. We denote by (x*, £„) the corresponding coordinates of the 
dual bundle 7r : E* — > M. One can also say that £ a is the fiber-wise linear local function i(e a ) on E* 
corresponding to the local section e a of E. We have also adapted local coordinates: 

(x\y a , X 3,y b ) inTE, inTE*, 
(^,y a , Pj ,^) inT*£, (aJ%e o ,Pi,V 6 ) in T*E* . 

It is well known (cf. [TB]) that the cotangent bundles T*E and T*E* are examples of double vector 
bundles: 

T*E* — > E , T*E — *-E*. 

TIE* T 7T£ 7T 

V Y Y Y 

E* > M E > M 

The concept of a double vector bundle is due to J. Pradines [331 [31], see also [2TJ [T5] ■ In particular, all 
arrows correspond to vector bundle structures and all pairs of vertical and horizontal arrows are vector 
bundle morphisms. The double vector bundles have been recently characterized |10j in a simple way 
as two vector bundle structures whose Euler vector fields commute. The above double vector bundles 
are canonically isomorphic with the isomorphism 

3v : T*E — > T* E* (2.1) 

being simultaneously an anti-symplectomorphism (cf. [El [12]). In adapted local coordinates, 3l T is 
given by 

%r(x\y a ,pj,TT b ) = (x l ,Tr a , -Pj,y b ). 

This means that we can identify Tib and as well as y b with (p b , so we will use local coordinates 
(x,y,p,£) on T*E and local coordinates (x,£,p,y) on T*E* in full agreement with the isomorphism 

For describing the dynamics of a Lagrangian or Hamiltonian systems it is necessary to give an 
additional ingredient, typically, a Lie algebroid structure or, more generally, the algebroid structure. 
In several papers, many authors have studied Lie algebroid structures as a unified geometric framework, 
general enough to account for different mechanical systems (defined on tangent bundles, on Lie algebras, 
quotients by Lie groups actions, etc.). 



It is well known that Lie algebroid structures on a vector bundle E correspond to linear Pois- 
son tensors on E* . A 2-contravariant tensor II on E* is called linear if the corresponding mapping 
II: T*E* — * TE* induced by contraction is a morphism of double vector bundles. This is the same 
as to say that the corresponding bracket of functions is closed on (fiber-wise) linear functions. The 
commutative diagram 

T*E* -^TE* , 




T*E 



composed with (2.1), describes a one-to-one correspondence between linear 2-contravariant tensors II 
on E* and homomorphisms of double vector bundles (cf. [H1II2]) covering the identity on E*: 
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T*E 



TE* 




TM 



In local coordinates, every e as above is of the form 

e(x\ y a ,p j7 Tr b ) = (x\ ir a , £ P 3 c (x)y c , £ ^abi^V^c + J2 a bi x )Pi) 



(2.2) 



(2.3) 



and it corresponds to the linear tensor 



a,6,c 



(2.4) 



i,6 



In |12j by algebroids we meant the morphisms (2.2 1 of double vector bundles covering the identity on 
E* , while Lie algebroids were those algebroids for which the tensor LT e is a Poisson tensor. We can 
consider the adjoint tensor LLj", i.e. the 2-contravariant tensor obtained from Il £ by transposition: 



a.b,c 



i.b 



and the corresponding adjoint algebroid structure e + . Algebroids e corresponding to skew-symmetric 
LI e (anti-symmetric brackets [-,-]e); i- e - such that e + = — e, we will call quasi-Lie algebroids. The 
relation to the canonical definition of Lie algebroid is given by the following theorem (cf. IT2"] L 

Theorem 2.1. j4n algebroid structure (E,s) can be equivalently defined as a bilinear bracket [•, -] £ on 
sections of t: E — > M, together with vector bundle morphisms p, a: E — > TM (left and right anchors), 
such that 

[fX, gY] e = fp(X)(g)Y - ga{Y){f)X + fg[X, Y} £ 

for f,g € C°°(M), X, Y G Sec(-E'). The bracket and anchors are related to the 2-contravariant tensor 
H £ by the formulae 

t([X,Y] s ) ={t(X),i(Y)}a e , 
n*(p(X)(f)) ={v(X),n*f} ne , 
n*(a(X)(f)) =U*f,L(X)}n s . 

The algebroid (E, e) is a quasi-Lie algebroid if and only if the tensor H E is skew-symmetric and it is a 
Lie algebroid if and only if the tensor Tl e is a Poisson tensor. 



The canonical example of a mapping e in the case of E = TM is given by e = Em 



the 



inverse to the Tulczyjew isomorphism olm ■ TT*M — > T*TM. In general, the algebroid structure map 
e is not an isomorphism and, consequently, its dual k — e* with respect to the right projection is a 
relation and not a mapping. Some authors (see [T3] for the case of Lie algebroids) have introduced the 
concept of prolongation of a Lie algebroid to maintain some of the original properties of the Tulzcyjew 
mapping (isomorphism and symplecticity) , but, for the purposes of this paper, it is not necessary to 
use this formalism. 
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The double vector bundle morphism (2.2 1 can be extended to the following algebroid analogue of 
the so called Tulczyjew triple 



T*E* 



TE* < 



T*E 



(2.5) 




The left-hand side is Hamiltonian, the right-hand side is Lagrangian, and the dynamics lives in the 
middle. 

We introduce now the dynamics through a Lagrangian L : E — > R which defines two smooth maps: 
the Legendre mapping: Al : E — > E* , Al = te* oeodL, which is covered by the Tulczyjew differential 
A L : E — > TE*, A L = e o dL: 




The lagrangian function defines the phase dynamics Tl = Al(E) C TE* which can be understood 
as an implicit differential equation on E* , solutions of which are 'phase trajectories' of the system 
(3 : R — > E* and satisfy t[3{t) e T L , where t/3 is the tangent prolongation of /?, t/3(i) = (/3(t), 0(t)). An 
analog of the Euler- Lagrange equations for curves 7 : R — » E is in turn: 

(£ L ) t(A L o 7 ) = A L o 7 . 



In local coordinates, Ll has the parametrization by (x l ,y a ) via in the form (cf. (2.3)) 
A L (x\y a ) = (x\ — (x, y),J2 pL(*)lf> E e c ab (x)y a ^- c (x, y) + ^ 4^)^, v)) 

" c a,c " i 

and the equation (El), for 7(f) = (x z (t), y a (t)), reads 



(2.6) 



(El) 



lit 



As one can see from (2.7 1, the solutions are automatically admissible curves in i.e. the velocity 
t(ro 7 )(t) is p( 7 (t)). 

With this framework it is possible to write in a unified point of view many equations of different 
mechanical systems that usually in the literature, appears as different ones (Classical Eulcr-Lagrange 
equations, Lagrange-Poincare equations after reduction by the action of a Lie group, Euler-Poincare 
equations, etc.). (See for applications of the theory in the case when (E,e) is a Lie algebroid). 



3. NONHOLONOMIC MECHANICS 

We can start from a general Lie algebroid (E, e) keeping in mind the standard case E = TM. A 
(linear) nonholonomic Lagrangian system is determined by a Lagrangian function L : E — > R and a 
vector subbundle D, rank!? = n — r, of the bundle E. We will accept subbundles over a submanifold, 
so let us denote Dm = t(D). By ijj : D <—* E let us denote the inclusion and by i* D : E? Dm — > D* the 
dual map. 

Because the solutions of the dynamics in E should be admissible curves, we need an admissibility 
condition ensuring that there are admissible curves through every point of D. The natural condition we 
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take is p(D) C TDm- By strong admissibility condition we will mean that D satisfies the integrabilty 
condition with respect to both algebroid structures: e and s + , i.e. 

p(D) C TD M and a(D) C TD M ■ (3.1) 

The Lagrangian function L : E — > K and the vector subbundle D define also the smooth map - the 
constrained Tulczyjew differential 

A£:D->T\D*, A% = Ti* D oeodL (3.2) 

covering the constrained Legendre map 

\1:D^D\ \® = i* D o \ L . (3.3) 

The diagram picture is the following 

T*£|D e - T(^ m ) T '° TD * (3.4) 

Af 



A 
dL 



Definition 3.1. The bundle D* we call the phase space of the nonholonomic system and the implicit 
differential equation L£ = A®(D), being a subset of TD*, we interpret as the phase dynamics of the 
nonholonomic system. A curve 7 : I D is a solution of the nonholonomic Euler-Lagrange equation 
if and only if A^ o 7 is an admissible curve in TD*, i.e. the nonholonomic Euler-Lagrange equation 
reads 

(nE L ): t(Af o 7 )=Af o 7 . (3.5) 



To find the explicit form of the nonholonomic Euler-Lagrange equation (3.5 1, consider local coordi- 
nates (x 1 ) — (x % , 1') on a open set U of M such that Dm is determined by the constraint x L = 0. A 
local basis {e a } a =i,...,n-r of sections of D we can extend to local sections of E and complete them to a 
local basis of sections {e a ,e a } of the vector bundle E. Then, in coordinates (x 1 ,y A ) — (x z , x L , y a 7 y a ) 
adapted to this bases, the local equations defining the constrained subbundle D as a vector subbundle 
of E over Dm are x L = 0, y a = 0. Note that admissibility of the constraint D means that p a (x l , 0) = 0. 

Taking local coordinates (x l ,y a ) on D, we may write then in : D <^-> E as iD(x l ,y a ) = (x l ,0,y a ,0) 
and i* D (x\0,£ a ,£ a ) = (x\£ a ), so 

^(^0,6^,0,6?) = (x*, . 

For the adapted local coordinates (x 1 ,x l ',y a ,y a ,Pj,p^,^b,S,i3) m T*E, the map e reduced to (T*E)m 
takes values in T(E^ D ^) (admissibility) and reads 

e{x\ 0, y a , 0,pj, Cb) - (x 4 , 0, £*, ^(i*, 0)/, 0, C^(^, 0)/£ c + ^(s*, 0>j) . 

Therefore 

TtJjoe^.Cj/'.O^j,^) = (x l ,^P^x\0)y d ,e^ b (x\0)y d ^c + ^(x\0)pj) (3.6) 



and 



Af(<z/ a ) 

dL dL 

= Ti* D (e(x\0,y a ,0 7 ^j(x\0,y a ,0),^(x\0,y a ,0)) 

fit fit fit 

= (x\ ^(x 1 , 0, y\ 0), p> A [x\ 0)y d , Q C db {x\ 0)/|^(x\ 0, y Q , 0) + a({x\ 0)^**, 0, y Q , 0)) 



Therefore, locally, the nonholonomic Euler-Lagrange equations read 

dt 

_ I dL^/^i ^ „,a I pC f^i n\„,d 9L / „i 

dt 



x'=0, y a = 0, d d f = p 3 d (x i ,Q)y d , (3.7) 
i(§^(x\0, y\ 0)) = t c db {x\ Q)y d ^b{x\ 0, y\ 0) + ^(x*, 0) 0, y Q , 0)) . (3.8) 
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In the case of a Lie algebroid a l A = p\, and if the subbundle D is over the whole base manifold 
M, Dm = M, the previous equations are precisely the nonholonomic equations obtained in [7] (see 
Equations 3.8). 

4. The nonholonomic reduction of algebroids 



It has been recognized a long time ago [22 35 that nonholonomic constraints may lead to certain 
nonholonomic brackets that do not satisfy the Jacobi identity. We will show that linear nonholonomic 
systems of mechanical type on general algebroids are again systems on general algebroids. But even if 
we start with a Lie algebroid, the new algebroid is, in general, no longer a Lie algebroid but certain 
quasi-Lie algebroid associated with a linear bi-vector field. This observation, made already in [S] 
(see also [HI), P u * s new n S n * *° the role of quasi-Poisson brackets, i.e. the brackets represented by 
arbitrary bi-vector fields and not satisfying, in general, the Jacobi identity (see [HI EH])- This shows, 
on the other hand, that developing Mechanics on general algebroids makes sense, as the reduction to 
a nonholonomic constraint will move us, in general, from the Lie algebroid picture into a more general 
one. 

Let e be an algebroid structure on a vector bundle E over M associated with the tensor IT e with the 



local form (2.4 1. For a linear subbundle D in E over a submanifold Dm C M, satisfying the strong 



admissibility condition (3.1 1, consider a decomposition 

E\ Dm = D @ Dm D x (4.1) 

and the associated projection P : Em M — ► D. We can construct the morphism ep : T*D — > TD* of 
double vector bundles as follows. 

To have the corresponding expressions in local coordinates, consider local coordinates (x 1 ) — (x l , x L ) 
on a open set U of M such that Dm is determined by the constraint x L — 0. Local bases {e a } a =i,..., n -r 
and {e c ,} Q=n _ r +i n of sections of D and D- 1 , respectively, can be extended to a basis &a of local 
sections of E. Then, we get the coordinates (x 1 ,y A ) — (x l , x\ y a , y a ) adapted to these bases and 
the adapted coordinates (x 1 \y A ,pj , £g) = (x z ,x l ',y a ,y a ,pj,p 7 ,^b,£,i3) in T*E. We get automatically 
coordinates (x l ,y a ) in D and (x l ,y a ,y a ) in Em M , and the adapted coordinates {x 1 , y a ,Pj, £ a ) and 
(x\y A ,p J ,^ B ) in T*D and T*(E {Dm ), respectively. 

Let us consider the phase lift T*P : T*D — \>T* (E\d m ) which, as often happens which phase lifts, 
is not a map but only a relation. In our local coordinates, 

T*P(x\y a 7 p r ^ b )^(x\y A , Pj ^ b ,0). 

We have also the embedding i^ |D : E\n M — > E whose phase lift, restricted to T*E\ 

T*(^ |D J : T*£p->T*(£ p J 

is the relation 

T*(i ElDM )(x\0,y a ,0 : pj,t B ) = (x\y A ,p,^ B ). 
The composition of relations 

T*(i Eln ) _1 oT*P:T*D->rEm 

has the local form 

(T*(i E]DM r 1 o T*P) (x\ y\p v f 6 ) = {x\ 0, y\ 0, Pj , 6,0). 



D ■ 



But, as o L h (x % , 0) = (strong admissibility condition), the value of Ti* D ° e in (3.6 1 does not depend on 
p-y and the composition 

e P = Ti* D oeo T*(iE lDM T 1 ° T*P : T*L> -> TD* (4.2) 

is a well-defined map which, in local coordinates, reads 

e P (x\ y a , Pj ^ b ) = (x\t a ,p> d (x\ 0)y d , G c db (x\ 0)y% + a k b (x\0)p k ) . (4.3) 

This is of course an algebroid structure on the bundle D. This algebroid structure can be described in 
a more straightforward way as follows. 
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The decomposition (4.1l gives the dual decomposition 



E \Dm = D * ®°m (D X T (4.4) 
and the corresponding projection P* : E'^ Dm — > D* . This projection does not give us a canonical 
projection of vectors and, more generally, contravariant tensors on E at points of E? Dm , unless they 
are tangent to B^. But the tensor LI £ with the local form (2.4 1 is a sum of tensor product with at 



least one part in the product being vertical. Moreover, due to strong admissibility condition for D, 
there is a unique decomposition Tl e = H £D + H £d± , where II eD is a linear tensor tangent to D* and the 
vertical parts of II e _,_ are tangent to (D- 1 )*, so they will be killed by the projection, independently 
how the other part of the product is. This gives a well-defined projection T(P*)(n e ) = H ep) when we 
restrict to the points of D* . 

To put it differently, we can take any smooth projection 7 from a neighbourhood of E? Dm onto 

E \Dm- Then ' 

T(P*)(II E ) :-T(FoT)(n £ ) (4.5) 
is a linear 2-contravariant tensor on D which does not depend on the choice of T. This tensor defines 
on D* the nonholonomic bracket {•, -} £p associated with the projection P. 

On the level of the algcbroid bracket [•, -] s this procedure is the following. The strong admissibility 
condition for D implies that the bracket [X, Y] e of sections of I? is a section of E\ Dm . Projecting this 
section to D along P gives us a bracket 

[X,Y] ep =P[X,Y] e (4.6) 

on sections of D - the nonholonomic restriction of [•, •] along P. This is an algebroid bracket with the 
original anchors. 

Of course, if the constraint bundle I? is a subalgebroid of E, then it is a holonomic constraint and 
the projection P plays no role: the tensor II e is tangent to D and we just take the restriction. It 
is however clear that, in the true nonholonomic case, the nonholonomic bracket need not satisfy the 
Jacobi identity, even when the original bracket does. In local coordinates, 

n eP = ® % + E eitf, °)% ® ^ - E o)^ ® d ta . 

a,b,c j.b a.j 

and the corresponding nonholonomic bracket reads: 

{x k ,x j } ep = 0. 

In the case of a Lie algebroid, we have & c ab = — C£ a , a J A = f? A , and the above bracket corresponds to 
the one introduced by [35] (see also [7j and references therein) for the subbundles over the total base 
M. 

A natural question arises here: what Lie algebroid structures we can obtain as nonholonomic re- 
strictions of Lie algebroid brackets? Of course, as the brackets must be skew-symmetric automatically, 
the algebroids must be necessarily quasi-Lie algebroids. We will call them nonholonomic quasi-Lie 
algebroids. In other words: which linear bi-vector fields on a vector bundle D — > M can be obtained 
by projections of Poisson tensors from a bigger vector bundle E — > M. For nonholonomic quasi-Lie 
algebras the answer is simple. 

Theorem 4.1. Any quasi-Lie algebroid D with the trivial anchor is a nonholonomic quasi-Lie alge- 
broid. 

Proof. As the anchor is trivial, we can clearly reduce to finite-dimensional quasi-Lie algebras. Take a 
basis ej of D and the corresponding linear coordinates y 1 , i — 1, . . . , n. The algebroid bracket in D is 
then determined by the structure constants Cy, 
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satisfying = — c^. Let us consider a new algebra E with the basis ej,/j, z = 1, . . . ,n and the 
bracket [-,-]e f° r which fj are central elements and 

[eue^E = c^f k . 

Since the algebra is 2-step nilpotent, it satisfies the Jacobi identity and it is a Lie algebra. We can 
them embed D in E by putting <p(ei) = (e, + /,) and take the complementary subspace D 1 - as spanned 
by elements a — /». The projection P : E ^ D is therefore given by 

^(ei) - e 4 + P(/i) - e t + /, . 

The nonholonomic bracket on the embedded submanifold D is therefore 

[et + fi,ej + fj] = P{c%f k ) = 4(e fe + f k ) , 

thus the original bracket in D. □ 



For a general quasi-Lic algcbroid the situation is much more complicated and we do not know a full 
characterization of nonholonomic quasi-Lie algebroids. Note however that not all quasi-Lie algcbroids 
are nonholonomic reductions of Lie algebroids. 

Example 4.2. Take a quasi-Lie algebroid structure e on D = TM with the anchor p = idTM- The 
corresponding tensor on T* M has therefore the local form 

n e = ^Cij(x)p k d Pi A d Pj + d Pi A d x i . 

Suppose that the corresponding bracket is a nonholonomic reduction of a Lie algcbroid bracket [-,-]e 
of a bigger vector bundle E along a projection P, 

[X,Y] S = P([X,Y] E ). 

Since the anchor p of D is the restriction of the anchor ps to D, the Lie algebroid E satisfies pe(E\m) = 
TM, thus E\m = D ®m K, where K is the kernel of Pe over M. For sections X, Y of D we have 
therefore pe([X,Y]e) = p{[X,Y]), thus, as the anchor map is for Lie algebroids a homomorphism of 
brackets, 

p([X,Y}) =p E ([X,Y] E ) = [p E (X),p E (Y)] vf = [p(X),p(Y)] vf , 

where [•, -] v f is the bracket of vector fields. We get that the anchor p maps the algcbroid bracket into 
the bracket of vector fields which is not the case for generally non-zero structure functions cjj- (x) . 

5. Nonholonomic systems of mechanical type 
Let us consider now on our general algebroid (E,e) a Lagrangian L : E — > M of mechanical type, i.e. 

L(e)= 1 -G(e,e)-V(r(e)), (5.1) 

where S : E xjj E — > R is a bundle metric on E and V : Ad — > R is the potential function. Let 
us consider also a a vector subbundlc D of E over satisfying the strong admissibility condition. 
Having the metric S in E, we have the natural decomposition E\ Dm — D ®d m w i tn being 
the orthogonal complement of D with respect to the metric £5, accompanied with the corresponding 
(this time - orthogonal) projection P : E\ Dm — > D. The fundamental observation in this case is the 
following. 

Theorem 5.1. The nonholonomic Tulczyjew differential A£ coincides with the Tulczyjew differential 
A; associated with the restricted Lagrangian I = L\ D : D — > M and the nonholonomic algebroid structure 
Ep on D. In other words, the nonholonomic Euler -Lagrange equations associated with L and the 
nonholonomic constraint D in E coincide with the non- constrained Euler- Lagrange equations for the 
Lagrangian function I on the algebroid (D,ep). 
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Proof. Let us consider local coordinates (x 1 ) = (x l : x L ) on a open set U of M such that Dm is 
determined by the constraint x L — 0. Take a local basis of orthonormal sections {e a , e a } of E adapted 
to the orthogonal decomposition E\ Dm = D ®\d m D , i.e, span{e a } = D and span{e Q } = D . We 
get the induced coordinates (x 1 ,y A ) in E in which the Lagrangian takes the form 



2 , 



Note that the 'mass' is 1 in these coordinates. On E\ Dm we have the induced coordinates (x z , y a , y a ) , 
so D is locally determined by the constraints x L — 0,y a — and the restricted Lagrangian is: 



l(x\y a ) = li2(y a ) 2 ~V(x\0). 



2 

a=l 

In our coordinates, the projector P reads 

P(x i ,y a ,y a ) = (x i ,y a ) 

and we deduce that 

s P (x\ y a , Pj , &) = (x\ C , (P d (x\ 0)y d , G c db (x\ 0)y% + a k (x\ 0) Pk ) 

so that 



e P o dl(x\ y a ) = (x\ y a ,p^x\0)y d 1 ^ Q c db {x\ 0)y d y c - a 



b dx k 



We get therefore the Euler-Lagrange equations for the Lagrangian I on D of the form 

which are exactly the nonholonomic equations ( |3.7| and ( |3.8| , if we take to account that 

riJ 

— (x\0,y a ,0) = SCy a 

and 

^(x\0) = 0. 

□ 

The above theorem makes it clear that, at least for Lagrangians of mechanical type, the nonholo- 
nomic Euler-Lagrange equations are just Euler-Lagrange equations, but on the nonholonomic restric- 
tion of the algebroid. In any case, passing to a liner nonholonomic constraint does not move us out 
of the Mechanics on algebroids. Our theory is therefore complete with respect to passing to the non- 
holonomic case, what was the main problem in understanding the nonholonomic constraint in the Lie 
algebroid (thus canonical) case. 

Having interpreted the nonholonomic Euler-Lagrange equations as just Euler-Lagrange equations, 
but on a reduced algebroid, we can use the generalization of Noether Theorem formulated for general 
algebroids in |Hj to obtain a nonholonomic Noether Theorem, at least for Lagrangians of mechanical 
type. 

The generalized Noether Theorem for a general algebroid (E, e) is based on the concept of the 
complete lift dj(K) of tensor fields K being sections of the tensor products E® k to the corresponding 
contravariant tensor fields on E. This is a natural generalization of the standard concept of the tangent 
lift dj which lifts contravariant tensors on a manifold M to the corresponding tensor fields on TM, 
(cf. HUI37]). 

For a vector bundle r : E — > M, let (8) fe (r) be the space of sections of the tensor-product bundle 
E® k over M. With any tensor field K G ® fc (r) we can associate the linear function t(K) on the dual 
bundle (E® k )* = (E*)® k and the vertical lift v T {K) € ^(^e) In local coordinates, 

i(f ai - ak (x)e ai (8 • • • (8 e a J = /«*~ * (x)U-a k 



NONHOLONOMIC CONSTRAINTS: A NEW VIEWPOINT 



11 



and 

A particular case of the vertical lift is the lift \ij(K) of a contravariant tensor field K on M into a 
contravariant tensor field on TM. It is well known (see [37J[TT]) that in the case of E = TM we have 
also the tangent lift dj: ® (tm) - * ®(ttm) which is a vj-derivation. What has been done for the 
tangent bundle, can be repeated in the case of an arbitrary algebroid (E,e). Note first that we can 
extend e naturally to mappings (cf. [T21 H2] ) 

e® r : ® r E T*E -^T® r M E*, r > 0. 

Theorem 5.2. [13l [12] Let (E,e) be an algebroid. For K £ ® k (r), k > 0, the equality 

L{&^{K))=& T {L{K))oe® k (5.2) 

defines the tensor field dj(K) € ® (te) which is linear and the mapping 

d^: ® fe ® fe (r) — ► © fc ® k (t e ) 

is a Vj -derivation of degree 0. In local coordinates, the lifts of functions on M and sections of E read: 

dHf(x)) = y- p i( x )^l( x ), (5.3) 

d|(r(x)e a ) = r(x)<7i(x)d xi +(y a pi(x)^(x) + ei b (x)y a f\x)]d y a. (5.4) 

Conversely, if D : (Bk ® (t) — ► ©& ® fc (rg) is a ^-derivation of degree smc/i f/iai D(K) is linear 
for each K G ® (t), £/ien i/iere is an algebroid structure s on r: E — > M suc/i i/ia£ Z? = dy. TTiis 
algebroid structure is a Lie algebroid if and only if 

d\{[X,Y] e ) = [d\{X)^{Y)] vf 

for all X,Y E ® 1 (r)- 

Suppose now that we are dealing with a Lagrangian L of mechanical type and a nonholonomic 
constraint subbundle D of (E,e) satisfying the strong admissibility assumption, so that we have also 
the nonoholomic reduced algebroid (D,ep). We will say that a pair (X, f) consisting of a section X of 
D and a function / on Dm is a symmetry of the nonholonomic problem associated with the Lagrangian 
L, if 

d?(X)(l)=d¥V), (5.5) 
where I = L\ D is the restriction of the Lagrangian to the constraint. 

Theorem 5.3. [Nonholonomic Noether Theorem] The following are equivalent: 

(a) the pair (X, f) is a symmetry of the nonholonomic problem associated with the Lagrangian I; 

(b) the function lx — f T D* on D* is a constant of the motion of the nonholonomic phase 
dynamics on D* , i.e. dy(tx — / ° To*) vanishes on r£; 

(c) the function 

{iX-fo T D *) ° Af 

is a constant of the motion for the nonholonomic Euler-Lagrange equation. 
Proof. The proof is completely analogous to that for Theorem 4 in [9] . □ 



For a section X of D, not verifying necessarily condition (5.5), we obtain 

d?(J0(Qo 7 (i) = ^( 4 *oA£o 7 (t)) (5.6) 

for any solution 7 : R — > D of the nonholonomic Euler-Lagrange equation. This last equation can be 
interpreted as a general version of the nonholonomic momentum equation studied for several authors 
(see [7J [2] and references therein) for nonholonomic systems with symmetry. 
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6. Examples 

Example 6.1. The Chaplygin sleigh. 

As an example of nonholonomic system on a Lie algebra, we study the Chaplygin sleigh which 
describes a rigid body sliding on a plane. The body is supported in three points, two of which slides 
freely without friction while the third point is a knife edge. This imposes the constraint of no motion 
orthogonal to this edge (see [5} 128]). 

The configuration space before reduction is the Lie group G = SE(2) of the Euclidean motions of 
the 2-dimensional plane M 2 . We will need in the sequel to fix some notation about the Lie algebra 
se(2). First of all its elements are matrices of the form 




and a basis of the Lie algebra se(2) = K 3 is given by 
Ex 





1 











(: 


0, 


r 

E 2 = 





1 


), *-( 




0/ 


V o 





oy 






[e 3 ,e ± ] 


= E 2 , [E 2 , E 3 ] 




Ei, 


[E U E 2 ] = 




We have that 
An element £ € se(2) is of the form 

£ = Vi Ex + v 2 E 2 + u> E 3 . 

The Chaplygin system is described by the kinetic Lagrangian function 

L : se(2) — ► R 

(vi, 1)2,0}) ' — ► | [( J + m(a 2 + b 2 ))uo 2 + mv\ + mv 2 — 2bmujv 1 — 2omwv 2 ] 

where m and J denotes the mass and moment of inertia of the sleigh relative to the contact point 
and (a, b) represents the position of the center of mass with respect to the body frame determined 
placing the origin at the contact point and the first coordinate axis in the direction of the knife axis. 
Additionally, the system is subjected to the nonholonomic constraint determined by the linear subspace 
of se(2): 

D = {{v u v 2 ,u)) G se(2) \ v 2 = Q}. 
Instead of {Ex, E 2 , E 3 } we take the basis of se(2): 

{ex — E 3 , e 2 = Ex, e 3 = —maE 3 — mabEx + ( J + ma 2 )E 2 } 

which is a basis adapted to the decomposition D © D ; D = span {ex, e 2 } and D — span {e 3 }. 
In the induced coordinates (jj 1 ,y 2 ) on D the restricted lagrangian is 

l(y\y 2 )= \ [(J + m(a 2 + b 2 ))(y 1 ) 2 + m (y 2 ) 2 ~2bmy 1 y 2 ] , 



and moreover, 



ma mab 
[ex,e 2 \ Ep = — + T~i 2 e2 ' 



Therefore, C} 2 = and C? 2 = -r^h. 

5 12 J+7na 2 12 J+ma 2 

Then, 

£p o dl(y , y )) — ((J + m(a 2 + fo 2 ))?/ 1 — bmy 2 ,my 2 — bmy 1 , —may 1 y 2 ,ma(y 1 ) 2 ) 

In consequence, the equations of motion are: 

(J + m(a 2 + b 2 ))^ 1 — bmy 2 = —may 1 y 2 
my 2 — bmy 1 — mafy 1 ) 2 



NONHOLONOMIC CONSTRAINTS: A NEW VIEWPOINT 



13 



ma 



ry 1 W-y 2 ) 



7 i 9 ' 

J + ma z 

y 2 = -T^^ymJ+mia' + b^-by 2 ) 



T I 2 ' 

J + mar 

Example 6.2. TTie snakeboard 

As a mechanical system the snakeboard has as configuration space Q — SE(2) x T 2 with coordinates 
(x,y,6,ip,<t>) (see [HEUrf]). 




The nonholonomic dynamics is described by 
• The Lagrangian 



= -m(i 2 + y 2 )+ (J + 2J 1 )0 2 + -J o (0 + ^) 2 + Ji i2 



2 x " ' 2 V y 2 

where m is the total mass of the board, J > is the moment of inertia of the board, Jo > 
is the moment of inertia of the rotor of the snakeboard mounted on the body's center of 
mass and J\ > is the moment of inertia of each wheel axles. The distance between the 
center of the board and the wheels is denoted by r. For simplicity (see [T7]), we assume that 
J + J + 2Ji = mr 2 . 

• The nonholonomic constraints induced by the non sliding condition in the sideways di- 
rection of the wheels: 

—x sin((9 + (/))+ y cos((9 + </>) — rO cos = 

—x sin((? — 4>) + y cos((9 — (p) + rO cos = 0. 

Observe that the Lagrangian is induced by the riemannian metric S on Q, 

9 = mdx 2 + mdy 2 + mr 2 dd 2 + J d9 eg) rp + J di[> <£> dO + J Q dijj 2 + 2J 1 dcf) 2 . 

The constraint subbundle td ■ D i — ► Q is 

_ i d d d d d 

D = span <^ ei = jrj,e 2 = ^r;,e 3 = a— + 6— + c— 
L a^y a0 ox ay 09 

where 

a = — r(cos 0cos((9 — 0) + cos 0cos((9 + <j))) — — 2r cos 2 0cos 
6 = — r(cos0sin(6' — 0) + cos0sin((9 + 0)) = — 2rcos 2 0sin(9 
c = sin(20). 

The orthogonal complement of D is spanned by 

{d d d d d 

e 4 = -b— + a—,e 5 — (cj — cmr 2 ) — + am— - am— 
ax ay ox o9 dip 

In the induced coordinates (x, y, 9, tp, 0, y 1 , y 2 , y ) on D the restricted lagrangian is 

l((x, y, 9, ip, 0, y\y 2 , y 3 ) = 2mr 2 cos 2 <p(y 3 ) 2 + J Q cy l y* + l ~My 1 ) 2 + My 2 ) 2 . 
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where now the nonholonomic constraints are rewritten as: y 4 = and y 5 — 0. After some straightfor- 
ward computations we deduce that 

[ei,e 2 ] £p = 0, 
[ei,e 3 ] £p = 0, 

2mr cos (f> (mr + cos 2cf>) tan 4> 

[e2,e 3 \ Ep = . 2 ei — e 3 . 

mr z — Jo sin <p mr z — Jo sin M 

Thus, the unique non vanishing structure functions are: 



2mr 2 cos 2 4> ra3 „ 3 (mr 2 + cos 2<p) tan <j> 



pi _ pi _ y P 3 _ P 3 _ 

^23 ~ — U 32 ~ 2 7 . 2 / ' ^-23 — — ^ 3 2 — ~ 2 7 • 2 

mr z — Jo sin q~> mr z — Jq sin ■ 



Therefore, the equation of motion of the snakeboard are 

x = ay 3 , y = by 3 , = cy 3 , ip = y 1 , <f> = y 2 , 

|(2/ 1 + cy 3 ) = 0, i(y 2 ) = 0, 
jt(4mr 2 cos 2 0y 3 + Joey 1 ) = 2J cos(2q>)y 1 y 2 - 2mr 2 sm(2(j))y 2 y 3 
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